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Abstract 



The beta-encoder was recently proposed as a quantization scheme for analog-to- 
digital conversion; in contrast to classical binary quantization, in which each analog 
sample x G [—1, 1] is mapped to the first N bits of its base-2 expansion, beta-encoders 
replace each sample x with its expansion in a base (3 between 1 < /3 < 2. This expan- 
sion is non-unique when 1 < (3 < 2, and the beta-encoder exploits this redundancy to 
correct inevitable errors made by the quantizer component of its circuit design. 

i-^h The multiplier clement of the beta-encoder will also be imprecise; effectively, the 

true value /3 at any time can only be specified to within an interval [Pi ow , Phigh] ■ 
This problem was addressed by the golden ratio encoder, a close relative of the beta- 

i— i encoder that does not require a precise multiplier. However, the golden ratio encoder 

is susceptible to integrator leak in the delay elements of its hardware design, and this 
has the same effect of changing /3 to an unknown value. 

In this paper, we present a method whereby exponentially precise approximations 
to the value of f3 in both golden ratio encoders and beta encoders can be recovered 
amidst imprecise circuit components from the truncated beta-expansions of a "test" 
number Xt es t G [ — 1, X] , and its negative counterpart, —x tes t- That is, beta-encoders 

\& and golden ratio encoders are robust with respect to unavoidable analog component 

< ^ imperfections that change the base (3 needed for reconstruction. 

00 

Keywords: A/D conversion, beta-encoders, beta-expansions, golden ratio encoder, 
quantization, redundancy, sampling, sigma-dclta 

c3 1 Introduction 



Beta-encoders with error correction were recently introduced in [3] as quantization algo- 
rithms in analog-to-digital (A/D) conversion that are simultaneously robust with respect 
to quantizer imperfections (like Sigma-Delta schemes) and efficient in their bit-budget use 
(like Pulse Code Modulation schemes). For a detailed discussion comparing Sigma Delta 
and Pulse Code Modulation (PCM) quantization schemes, we refer the reader to [6] and 
[3]. Recall that after the sampling step in analog-to-digital conversion of a bandlimited 
function f(t) with L°° norm ||y||oo — 1? *^n N-bit PCM quantizer simply replaces each 
sampled function value x = f(t n ) with the first N bits of its truncated binary expansion. 
Beta-encoders are similar to PCM; they replace each sampled function value x with a 

*R. Ward is a graduate student in the Program in Applied and Computational Mathematics, Prince- 
ton University, Princeton, New Jersey 08544 (e-mail: rward@princeton.edu). This research project was 
supported by an NSF graduate research fellowship. 
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truncated series expansion in a base 0, where 1 < < 2, and with binary coefficients. 
Clearly, if = 2, then this algorithm coincides with PCM. However, whereas the binary 
expansion of almost every real number is unique, for every € (1,2), there exist a con- 
tinuum of different expansions of almost every x in (0,1] (see |8]). It is precisely this 
redundancy that gives beta-encoders the freedom to correct errors caused by imprecise 
quantizers shared by Sigma Delta schemes. Whereas in Sigma Delta, a higher degree of 
robustness is achieved via finer sampling, beta-encoders are made more robust by choosing 
a smaller value of as the base for expansion. 

Although beta-encoders as discussed in [3] are robust with respect to quantizer imperfec- 
tions, these encoders are not as robust with respect to imprecisions in other components 
of their circuit implementation. Typically, beta-encoders require a multiplier in which real 
numbers are multiplied by 0. Like all analog circuit components, this multiplier will be 
imprecise; that is, although a known value 0$ may be set in the circuit implementation 
of the encoder, thermal fluctuations and other physical limitations will have the effect of 
changing the true multiplier to an unknown value (3 £ [Plow > Phigh] within an interval of 
the pre-set value 0q. The true value (3 will vary from device to device, and will also change 
slowly in time within a single device. This variability, left unaccounted for, disqualifies the 
beta-encoder as a viable quantization method since the value of (3 must be known with 
exponential precision in order to reconstruct a good approximation to the original signal 
from the recovered bit streams. 

We overcome this potential limitation of the beta-encoder by introducing a method for 
recovering (3 from the encoded bitstreams of a real number x G [—1,1] and its negative, 
— x. Our method incorporates the techniques used in [2], but our analysis is simplified 
using a transversality condition, as defined in [1], for power series with coefficients in 
{—1,0, 1}. As the value of can fluctuate within an interval [0iow,0high] over time, our 
recovery technique can be repeated at regular intervals during quantization (e.g., after the 
quantization of every 10 samples). 

The golden ratio encoder (GRE) was proposed in [5 a as a quantizer that shares the same 
robustness and exponential rate-distortion properties as beta-encoders, but that does not 
require an explicit multiplier in its circuit implementation. GRE functions like a beta- 
encoder in that it produces beta-expansions of real numbers; however, in GRE, is fixed 
at the value of the golden ratio, = (f> = 1+ 2 V ^ . The relation <j) 2 = <p + 1 characterizing 
the golden ratio permits elimination of the multiplier from the encoding algorithm. Even 
though GRE does not require a precise multiplier, component imperfections such as in- 
tegrator leakage in the implementation of GRE may still cause the true value of to be 
slightly larger than (f>; in practice it is reasonable to assume £ [4>, 1.1(f)]. Our method for 
recovering in general beta-encoders can be easily extended to recovering in the golden 
ratio encoder. 

The paper is organized as follows: 

1. In sections I (a) and I (b), we review relevant background on beta-encoders and 
golden ratio encoders, respectively. 

2. In section II, we introduce a more realistic model of GRE that takes into account 
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the effects of integrator leak in the delay elements of the circuit. We show that the 
output of this revised model still correspond to truncated beta-expansions of the 
input, but in an unknown base (3 that differs from the pre-set value. 

3. Section III describes a way to recover the unknown value of f3 up to arbitrary pre- 
cision using the bit streams of a "test" number x G [—1, 1], and —x. The recovery 
scheme reduces to finding the root of a polynomial with coefficients in {—1,0, 1}. 

4. Section IV extends the recovery procedure of the previous section to the setting of 
beta-encoders having leakage in the (single) delay element of their implementation. 
The analysis here is completely analogous to that of Section III. 

1.1 A brief review of beta-encoders 

In this section, we summarize certain properties of beta-encoders with error correction, 
from the perspective of encoders which produce beta expansions with coefficients in 
{— 1, 1}. For more details on beta-encoders, we refer the reader to [3]. 

We start from the observation that given f3 G (1,2], every real number x G [— ^ry, pzj] 
admits a sequence (bj)j^N, with bj £ {—1,1}, such that 

oo 

Under the transformation bj = -^j— , is equivalent to the observation that every real 
number y G [0, ^-j-] admits a beta-expansion in base /3 G (1, 2], with coefficients bj G {0, 1}. 
Accordingly, all of the results that follow have straightforward analogs in terms of {0, 1}- 
beta-expansions; see [3] for more details. 

One way to compute a sequence (bj)j^N that satisfies ([I]) is to run the iteration 

f3x 

Q{ui) 

P(uj - bj) 

Q(u j+1 ) (2) 
where the quantizer Q(u) is simply the sign- function, 

For j3 = 2, the expansion ([I]) is unique for almost every x G [—1, 1]; however, for (3 G (1, 2), 
there exist uncountably many expansions of the type ([I]) for any x G [—1,1] (see [8]). 
Because of this redundancy in representation, beta encoders are robust with respect to 
quantization error, while PCM schemes are not. We now explain in more detail what 
we mean by quantization error. The quantizer Q(u) in Q is an idealized quantizer; in 
practice, one has to deal with quantizers that only approximate this ideal behavior. A 



Ul = 

h = 

for j > 1 : Uj+i = 

bj+i = 
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more realistic model is obtained by replacing Q(u) in ([3]) with a "flaky" version Q u (u), 
for which we know only that 

u < —v 

u>v (4) 
— V <u<v. 

In practice, v is a quantity that is not known exactly, but over the magnitude of which we 
have some control, e.g. \v\ < e for a known e. This value e is called the tolerance of the 
quantizer. We shall call a quantization scheme robust with respect to quantization error if, 
for some e > 0, the worst approximation error produced by the quantization scheme can be 
made arbitrarily small by allowing a sufficiently large bit budget, even if the quantizer used 
in its implementation is imprecise to within a tolerance e. According to this definition, 
the naive {—1, l}-binary expansion is not robust. More specifically, suppose that a flaky 
quantizer Q u (u) is used in ^ to compute the base-2 expansion of a number x G [—1,1] 
which is sufficiently small that \2x\ < v. Since 2x is within the "flaky" zone for Q u (u), if 
b\ = Q u {2x) is assigned incorrectly; i.e., if b\ differs from the sign of x, then no matter 
how the remaining bits are assigned, the difference between x and the number represented 
by the computed bits will be at least \x\. This is not the case if 1 < (3 < 2, as shown by 
the following theorem whose proof can be found in [3] : 

Theorem 1.1. Let e > and x G [—1,1]. Suppose that in the beta- encoding of x, the 
procedure ([2| is followed, but the quantizer Q v {u) is used instead of the ideal Q(u) at each 
occurence, with v satisfying v < e. Denote by (bj)j^jsf the bit sequence produced by applying 
this encoding to the number x. If (3 satisfies 

then 

N 

\x-^2bjP- j \ <C(3~ N (5) 

with C = e + 1 . 

For a given tolerance e > 0, running the recursion ^ with a quantizer Q u {u) of 
tolerance e and a value of (3 in (1, produces bitstreams (bj) corresponding to a beta- 
expansion of the input x in base (3; however, the precise value of (3 must be known in order 
recover x from such a beta-expansion. As mentioned in the introduction and detailed in 
section II, component imperfections may cause the circuit to behave as if a different value 
of j3 is used, and this value will possibly be changing slowly over time within a known range, 
[Plow , Phigh] ■ In [2] , a method is proposed whereby an exponentially precise approximation 
7 to the value of 7 = at any given time can be encoded and transmitted to the 
decoder without actually physically measuring its value, via the encoded bitstreams of a 
real number x G [0, 1) and 1 — x. This decoding method can be repeated at regular time 
intervals during the quantization procedure, to account for the possible time- variance of 7. 
That an exponentially precise approximation 7 to 7 is sufficient to reconstruct subsequent 
samples f(t n ) with exponential precision is the content of the following theorem, which is 
essentially a restatement of Theorem 5 in [2] . 

Theorem 1.2 (Daubechies, Yilmaz). Consider x £ [0,1) and (bj)jzN G {0)1}; or x £ 
[— 1, 1] and (bj)j£N G {—1, 1}. Suppose 7 G (1/2, 1) is such that x = YlTLi ■ Suppose 



Q v {u) 
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7 is such that \j — 7I < C\j 



for some fixed C\ > 0. Then xjy : = Ylj=i bjrf satisfies 



x - x N \ < C27 



N 



(6) 



where C2 is a constant which depends only on 7 and C\. 

Although the approach proposed in [2] for estimating (3 from bitstreams overcomes 
potential approximation error caused by imprecise multipliers in the circuit implementa- 
tion of the beta-encoder, new robustness issues are nevertheless introduced. Typically, 
one cannot ensure that the reference level 1 in 1 — x is known with high precision. To 
circumvent this problem, the authors consider other heuristics whereby (3 is recovered 
from clever averaging of multiple pairs Xj and 1 — Xj. These heuristics do not require 
that the reference level 1 in 1 — x be precise; however, these approaches become quite 
complicated in and of themselves, and any sort of analytical analysis of their performance 
becomes quite difficult. In section III, we will present a similar approach for recovering 
(3 that does not require a precise reference level, yet still allows for exponentially precise 
approximations to (3. 

1.2 The golden ratio encoder 

As shown in the previous section, beta-encoders are robust with respect to imperfect 
quantizer elements, and their approximation error decays exponentially with the number 
of bits N. To attain this exponential precision, (3 must be measured with high precision, 
which is quite complicated in practice. These complications motivated the invention of 
the golden ratio encoder (GRE) of [5 , which has the same robustness and rate-distortion 
properties as beta-encoders, but uses an additional delay element in place of precise mul- 
tiplier in its implementation. More precisely, if one implements the recursion Q with 
(3 = 4'= 1+ 2 V ^ , then using the relation (ft 2 = (ft + 1, one obtains the recursion formula 
u n+ \ = u n _i + u n — (b n -i + (ftb n ). If the term b n -\ + (ftb n in this formula is removed, 
then the resulting recursion v n+ \ = v n + v n -i should look familiar; indeed, with initial 
conditions (vq,Vi) = (0,1), this recursion generates the Fibonacci numbers v n , and it is 
well-known that the sequence — > (ft as n — > 00. If 6 n _i + <ftb n is instead replaced by a 
single bit taking values in {—1,1}, then we are led to the following scheme: 



"0 = 



U\ = x 




(7) 




u + av < 
u + av > 



(8) 



along with their flaky analogs 




-1, 
1, 

— 1 or 1 



u + av < —v 

u + av > v 

—v < u + av < v 



0) 
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In [S], the authors consider the recursion formula ([7]) implemented with flaky {0,1}- 
quantizers of the form (u,v) = Q u (u + av — l) + 1 /2. For the simplicity of pre- 
sentation, we will consider only the { — 1, l}-quantizers Q, but many of our results extend 
straightforwardly to quantizers of the type Q% L (u,v). 

The following theorem was proved in [5|; it shows that as long as x and Q(u, v) are such 
that the state sequence u = {« n }^L remains bounded, a golden ratio encoder (correspond- 
ing to the recursion d7|) will produce a bitstream (bj) corresponding to a beta-expansion 
of x in base /3 = (f>, just as does the beta-encoder from which the GRE was derived. 

Theorem 1.3 (Daubechies, Giintiirk, Wang, Yilmaz). Consider the recursion ([7]). Sup- 
pose the 1-bit quantizer Q(u,v) that outputs bits (bj) in { — 1,1} is of the type Q v a {u,v) 
such that the state sequence u = {u n }^ =0 with uq = and u\ = x is bounded. Then 

N 

\ x ~Y,Kr n \ <r N+1 (io) 



n=l 



Here 4> is the golden ratio. 
Proof. Note that 

N N 



E & ™0 n = ^(«n-l + U n - U n+ i)(j) n 
n=l n=l 

N-l N N+l 

n=0 n=l n=2 

N-l 



= u^- 1 + + 1-0) E 

n=2 

+U 1 ((J)~ 1 + 0~ 2 ) + 0~ N (u N (l - Cj)) - UAT+l) 
= Uq^' 1 + Ul + (f)~~ N (u N (l - f) - U N+ ^j 

= x + 4>~ N (u N (l - (f>) - U N+ lj. 

The second to last equality uses the relation 4>~ 2 + (j)^ 1 — 1 = 0, and the last equality is 
obtained by setting no = and u\ = x. Since the state sequence u = {u n }^ =0 is bounded, 
it follows that x = YlnLi bn<P~ n ■ Thus, 



N 

E 

n=l n =N+l 



bn<t> n \ = | E hn( < 
--N+1 
(N+l) 



< 



1 - (j)" 1 
A-N+l 



(11) 



□ 
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(a) Beta-encoder (b) GRE-encoder 



Figure 1: A comparison of the block diagrams for (a) a beta-encoder with error correction, and 
a (b) golden ratio encoder. Here D represents a delay term, which holds input throughout one 
timestep. 



Figure 1 compares the block-diagrams of a beta-encoder with error correction, corre- 
sponding to the recursion formula ([2]), and a golden ratio encoder corresponding to the 
recursion ([7]). Although the implementation of GRE requires 2 more adders and one more 
delay element than the implementation of the beta-encoder, the multiplier element a in 
GRE does not have to be precise (see section 6), whereas imprecisions in the multiplier 
element (5 of the beta-encoder result in beta-expansions of the input x in a different base 

2 GRE: The Revised Scheme 

In modeling the golden ratio encoder by the system Q, we are assuming that the delay 
elements used in its circuit implementation are ideal. A more realistic model would take 
into account the effect of integrator leak, which is inevitable in any practical circuit im- 
plementation (see [9] for more details). After one clock time, the stored input in the first 
delay is reduced to Ai times its original value, while the stored input in the second delay 
is replaced by A2 times its original value. In virtually all circuits of interest, no more than 
10 percent of the stored input is leaked at each timestep; that is, we can safely assume 
that Ai and A2 are parameters in the interval [.9,1]. The precise values of these parameters 
may change in time; however, as virtually all practical A/D converters produce over 1000 
bits per second (and some can produce over 1 billion bits per second), we may safely as- 
sume that Ai and A2 are constant throughout the quantization of at least every 10 samples. 

Fixing an input value x £ [—1,1], we arrive at the following revised description of the 
GRE encoder: 

uo = 

Ul = X 

h = Q(0,x) 
for n > 1 : u n+ \ = \\u n + Ai A2U„,_i — b n 

K+l = Q(Un,Un-l) (12) 

Obviously, Ai = A2 = 1 corresponds to the original model Q. It is reasonable to assume 
in practice that (Ai, A2) G M := [.95, l] 2 , and in virtually all cases (Ai, A2) 6 V := [.9, l] 2 . 



We will show that the revised scheme (12) still produces beta-expansions of the input 
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x, but in a slightly different base 7 = (3 1 = ^^^x^^^ 2 ? which increases away from 



4>~ l as the parameters Ai and A2 decrease. Key in the proof of Theorem (10) was the use 
of the relation 1 — <j)~ l — (j)~ 2 = to reduce Y^n=l(. u n + u n+i — u n+ 2)<p~ n to the sum of the 
input x, and a remainder term that becomes arbitrarily small with increasing N. Accord- 
ingly, the relation 1 - A17 - A1A27 2 = gives J2n=i( u n + u n+i ~ u n +2h' n = x + R(N), 
where R(N) goes to as N goes to infinity. 



Theorem 2.1. Suppose the 1-bit quantizer Q in (j 1 2 p of type (|9|) is such that the state 

2A1A2 



sequence u = {u n }^ =0 with uq = and u± = x is bounded. Consider 7 = Al+ 2A^;it 4AlA2 
Then 



x 



En=l b nJ n \ < C^ N 



where C~ = 
Proof. 



N N 



^ &n7™ = ^(AiA 2 n n „i + XxUn - u„+i)7 r 



n=l ra=l 



JV-1 

= (AiA 27 + Ai - 7- 1 ) 7™^. 

n=2 

+AiA 2 u 7 + ^i(Ai7 + A1A27 2 ) 
+7 N (u N (Xi - 7 _1 ) - u N+1 ^j 

= x + 7 7V (n^(Ai- 7 ~ 1 )-^ +1 ). (13) 

The last equality is obtained by setting no = and u\ = x. Since the u n are bounded, it 
follows as in the proof of ([7]) that 

\x-J2bnl n \ < } • (14) 

1 — 7 

n=l ' 



□ 



Theorem (2.1) implies that if 7 is known, then the revised GRE scheme (12) still 
gives exponential approximations to the input signal x, provided that the u n are indeed 
bounded. The following theorem gives an explicit range for the parameters (u, a) which 
results in bounded sequences u n when the input x G [—1,1], independent of the values 
of the leakage parameters (Ai, A2) in the set V = [.9, l] 2 . This parameter range is only 
slightly more restrictive than that derived in [5] for the ideal GRE scheme ([7]); that is, 
the admissable parameter range for a and v is essentially robust with respect to leakage 
in the delay elements of the GRE circuit implementation. 



Theorem 2.2. Let x £ [-1, 1], and (Ai, A 2 ) G [.9, 1] . Suppose that the GRE scheme (12) 
is followed, and the quantizer Q^(u, v) is used, with v possibly varying at each occurrence, 
but always satisfying u < e for some fixed tolerance e < .5. // the parameter a takes values 
in the interval [1.198 + 1.198e, 2.281 — .952e], then the resulting state sequence {uj)j^M is 
bounded. 
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We leave the proof of Theorem (2.2) to the appendix. This theorem in some sense 
parallels Theorem (1.1) in that an admissable range for the "multiplication" parameter 
(a in GRE, (3 in beta-encoders) is specified for a given quantizer tolerance e; however, 
we stress that the specific value of (5 is needed in order to recover the input from the 
bitstream (&,) in beta-encoders. In contrast, a GRE encoder can be built with a multiplier 
a set at any value within the range [1.198 + 1.198e, 2.281 — .952e], and as long as this 
multiplier element has enough precision that the true value of a will not stray from an 
interval [ai ow , cthigh] C [1.198 + 1.198e, 2.281 — .952e], then the resulting bitstreams (bj) 
will always represent a series expansion of the input x in base 7 of Theorem (2.1 ), which 
does not depend on a. The base 7 does however depend on the leakage parameters Ai 
and A2, which are not known a priori to the decoder and can also vary in time from 
input to input: as discussed earlier, the only information available to the decoder a priori 
is that (Ai,A2) G [.9, l] 2 in virtually all cases of interest, and (Ai,A2) G [.95, l] 2 , or 
7 G jg^ -1 ] ~ [-618, .6505] in most cases of interest. In the next section, we show 

that the upper bound of .6505 is sufficiently small that the value of 7 can be recovered 
with exponential precision from the encoded bitstreams of a pair of real numbers (x, —x). 
In this sense, GRE encoders are robust with respect to leakage in the delay elements, 
imprecisions in the multiplier a, and quantization error. 



3 Determining 7 

3.1 Approximating 7 using an encoded bitstream for x = 



Recall that by Theorem (1.2), exponentially precise approximations 7 to the root 7 in 
Theorem (2.1 ) are sufficient in order to reconstruct subsequent input x n = f(t n ) G [—1, 1] 
whose bit streams are expansions in root 7 with exponential precision. In this section, we 
present a method to approximate 7 with such precision using the only information at our 
disposal at the decoding end of the quantization scheme: the encoded bitstreams of real 
numbers x G [—1, 1]. More precisely, we will be able to recover the value 7 = /3" 1 using 
only a single bitstream corresponding to a beta-expansion of the number 0. It is easy to 
adapt this method to slow variations of 7 in time, as one can repeat the following pro- 
cedure at regular time intervals during quantization, and update the value of 7 accordingly. 

The analysis that follows will rely on the following theorem by Peres and Solomyak [I]: 

Theorem 3.1 (Peres-Solomyak). [5-transversality] Consider the intervals I p = [0,p\. If 
p < .6491..., then for any g of the form 

00 

g{x) = l + ^bjaP, bje {-1,0,1} (15) 
i=i 

and any x G I p , there exists a 5 P > such that if g(x) < 5 P then g'(x) < —5 p . Furthermore, 
as p increases to .6491..., 5 P decreases to 0. 



Theorem (3.1) has the following straightforward corollary: 



Corollary 3.2. If g{x) is a polynomial or power series belonging to the class B given by 

00 

B = {±l + J2bjx j , '■ °j G {—1,0, 1}} (16) 
i=i 
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then g can have no more than one root on the interval (0, .6491]. Furthermore, if such a 
root exists, then this root must be simple. 



In pQ, Peres and Solomyak used Theorem (3.1) to show that the distribution v\ of 



the random series ^ ±A n is absolutely continuous for a.e. A G (1/2,1). The estimates 



in Theorem (3.1) are obtained by computing the smallest double zero of a larger class of 
power series B := {1 + Y^=i bn,x n '■ b n S [—1, 1]}. The specific upper bound p = .6491... 
for which 5-transversality holds on the interval I p is the tightest bound that can be reached 
using their method of proof, but the true upper bound cannot be much larger; in [7], a 



power series g(t) belonging to the class B in (16) is constructed which has a double zero 
at to « .68 (i.e., g(to) = and g'(to) = 0), and having a double root obviously contradicts 
5-transversality. 



We now show how to use Theorem (3.1) to recover 7 from a bitstream (bj) ! ^ =1 produced 
from the GRE scheme ( 12 ) corresponding to input x = 0. We assume that the parameters 
(a, v) of the quantizer Q^{u, v) used in this implementation are within the range provided 
by Theorem (2.2). Note that such a bitstream (bj)'jL 1 is not unique if v > 0, or if more 



than one pair (Ai, A2) correspond to the same value of 7. Nevertheless, Theorem ( |2.l[ ) im- 
plies that Yl'jLi ^j7 J = 0) so that 7 is a root of the power series F(t) = b\ + Yl'jLi bj+iP- 
Suppose we know that (3 > 1.54056, or that 7 < .6491. Since F(t) belongs to the class 



B of Corollary (3.2), 7 must necessarily be the smallest positive root of F(t). In reality 



one does not have access to the entire bitsream (bj)°^ =1 , but only a finite sequence {bj)^ =1 . 
It is natural to ask whether we can approximate the first root of F{t) by the first root 
of the polynomials P n {t) = b\ + J2j=i bj+it 1 . Since the P n still belong to the class B of 



Corollary (3.2), |P n (^)| has at most one zero on the interval [^-\.6491] » [.618, .6491]. 
The following theorem shows that, if it is known a priori that 7 < .6491 — e for some 
e > 0, then for n sufficiently large, |P n (i)| is guaranteed to have a root 7 n in [0, .6491], 
and 1 7 — 7„| decreases exponentially as n increases. 

Theorem 3.3. Suppose that for some e > 0, it is known that 7 < yhigh = 

.6491 - e. Let 

5 > be such that 5-transversality holds on the interval [O^high]- Let N be the smallest 
integer such that ^ N+l < (1 — j)e5. Then for n> N, 

(a) P n has a unique root 7„ in [0, .6491] 

(b) |7 - ln\ < Ci7 n , where C x = j^hfy 

Proof. Without loss of generality, we assume b\ = 1, and divide the proof into 2 cases: (1) 
Pn(l) < 0, and (2) P n (j) > 0. The proof is the same for b\ = —1, except that the cases 
are reversed. 

Case (1): In this case, many of the restrictions in the theorem are not necessary; in 
fact, the theorem holds here for all n > and e > 0. P n (t) has opposite signs at t = 
and t = 7, so P n must have at least one root 7™ in between. Moreover, this root is unique 



by Theorem (3.1). To prove part (b) of the theorem, observe that Theorem (3.1) implies 
that if P n (to) < S at some t G (0, .6491), then P' n (t) < -5 in the interval [t , .6491). 
In particular, P n (t) < 5 and P^if) < — 5 for t £ [7m 7]- By the Mean Value Theorem, 
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\P n (i)\ = \Pn{l) ~ P„(7«)| = |Pn(0ll7 " 7n| for some £ G [ 7n , 7 ], so that 

\Pn{l)\ 



17 - 7n 



< 

< 
< 



inf Ce[ 7 n,7] l P n(OI 
|P«(7)I 



5 

,,n+l 



(1-7)* 



v«+l 



The inequality |-Fn( 7 )| < ^_ follows from Theorem (2.1). 



Case (2): If P n {l) > 0, then Theorem (3.1) implies that the first positive root of P n , if 



it exists, must be greater than 7 . Whereas in Case (1), P n was guaranteed to have a root 
7n < 7 for all n > 0, in this case P n might not even have a root in (0, 1); for example, the 
first positive root of the polynomial P±{t) = 1 — t — t 2 + i 3 occurs at t = 1. However, if ra 
is sufficiently large, P n will have a root in [7, .6491]. Precisely, let n > N, where N is the 
smallest integer such that j N+1 < (1 — 7 )e<5. Then P n {l) < ~\~ — \~ < eS < 5. Since 



7 + e < .6491, Theorem (3.1 ) gives that P^(t) < -5 for t G [7, 7 + e], and 



Pn{l + e) < P„( 7 ) + e sup P' n {t) 

< e5 -eS 
= 0. 

We have shown that P n {j) > and P n (j + e) < 0, so P n is guaranteed to have a root 7n 
in the interval [7,7 + e], which is in fact the unique root of P n in [0, .6491]. Furthermore, 
the discrepency between 7 n and 7 becomes exponentially small as n > N increases, as 

\Pnh)\ 



U _ ~\ < < 

|7n 71 ^ inf^r.. , IP' mi - 



■ 

□ 

Remark 3.4. In p] it is shown that 5-transversality holds on [0, .63] with 5 = 5.63 = -07. 



This interval corresponds to e = .6491 — .63 = .0191 in Theorem (3.3). If 7 is known a 

7 

1-7 



priori to be less than .63, then ^— < 1 f i3 63 ~ 1.7027, and iV 0/ Theorem (3.3) satisfies 



M < log e<5-log 1.7027 <• 1fi 
iV - log .63 - iD - 

Of course, even if we know P n (t), n will be too large to solve for 7n analytically. 
However, if we can approximate 7n by 7n with a precision of 0(7"), then I7 — 7n | < 
1 7 — 7n\ + \ln — ln\ wm a ls° be of order 7 n , so that the estimates 7 n are still exponentially 
accurate approximations to the input signal. 

Indeed, any 7n G .6491 — e] satisfying |P n ( 7ri )| < <ft~ n provides such an exponen- 

tial approximation to 7 n . Since n > N, it follows that (j)~ n < 7™ < 7^ < — < 5, and 
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so Pn(t) < —5 on the interval between ^ n and j n . Thus, 

|7-7n| < 17 - 7n| + \ln ~ 7n| 

< c a n + 



Pn{ln)\ 



inf 



< C l7 n + 

< c a n + 
= c 2l n 



7" 



(17) 



where C 2 = Ci + i = 



3.2 Approximating 7 with beta-expansions of (— x,x) 

The method for approximating the value of 7 in the previous section requires a bitstream b 
corresponding to running the GRE recursion (12) with specific input x = 0. This assumes 



that the reference value can be measured with high precision, which is an impractical 
constraint. We can try to adapt the argument using bitstreams of an encoded pair (x, —x) 
as follows. Let b and c be bitstreams corresponding to x and —x, respectively. Define 
dj = bj + Cj. Put k = min{j £ N\dj + i ^ 0}, and consider the sequence d = (dj)JL l 
defined by dj = \dj + k- Since dj G {—2,0,2}, it follows that dj € {—1,0,1}, and by 



Theorem (2.1), we have that 



n=l 



^ 00 
L n=l 



Cn7 



n=l 



n=l 



so that 



iV 



(18) 



(19) 



where the constant C~ 



7 
1-7' 



Equation (19), along with the fact that the polynomials -Pjv(t) := d\ + Ylj=i dj+it-* 



are of the form (16), allows us to apply Theorem (3.3) to conclude that for N suffi- 



ciently large, the first positive root 77V of P/v becomes exponentially close to 7. However, 
note that the encoding of (x, —x) is not equivalent to the encoding of 0. The value of 
k = min{j £ N\dj + \ ^ 0} used to define the sequence d can be arbitrarily large. In 
fact, if an ideal quantizer Q^, is used, then the bitstreams b and c are uniquely defined by 
b = — c, so that d = 0. Thus, this method for recovering 7 actually requires the use of a 
flaky quantizer Q v a . To this end, one could intentially implement GRE with a quantizer 
which toggles close to, but not exactly at zero. One could alternatively send not only 
a single pair of bitstreams (6, c), but multiple pairs of bitstreams (b l ,c l ) corresponding 
to several pairs (x/, — x;), to increase the chance of having a pair that has b l ,c l / for 
relatively small j. 
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Figure 3 plots several instances of Ps,Pw, and P32, corresponding to 7 = .64375. The 
quantizer Q^(u,v) is used, with a = 2 and v = .3. These values of a and v generate 



bounded sequences (ti„) for all (Ai, A 2 ) G [.9, l] 2 by Theorem dgjb. 



As shown in Figure 2, numerical evidence suggests that 10 iterations of Newton's method 
starting from xq = 4>~ 1 ~ .618 will compute an approximation to 7/v, the first root of 
PjVj with the desired exponential precision. The figure plots 7 versus the error I7 — t/v|, 
where 7at is the approximation to 7jv obtained via a 10- step Newton Method, starting 
from xq = .618. More precisely, for each N, we ran 100 different trials, with x and 7 
picked randomly from the intervals [—1,1] and [.618, .7] respectively, and independently 
for each trial. The worst case approximation error of the 100 trials is plotted in each case. 
The quantizer used is Q^(it, v) with v = .3, and a picked randomly in the interval [1.7, 2], 



independently for each trial. Again, Theorem (2.2) shows that these values of a and v 
generate bounded sequences (u n ) for all (Ai, A2) G [.9, l] 2 . 

4 Bet a- expansions Revisited 

Even though our analysis of the previous section was motivated by leaky GRE encoders, 
it can be applied to general beta-encoders to recover the value of (5 at any time during 
quantization. From the last section, we have: 

Theorem 4.1. Let F(t) = Yl^jLo^jt 1 be a power series belonging to the class B = {±1 + 
Yj%.\ bj x3 bj G {-1,0, 1}}. Suppose that F has a root at 7 G [rflow,7high], where j high = 
.6491 — e for some e > 0. Let S > be such that 5-transversality holds on the interval 
[0,jhigh\- Let N be the smallest integer such that ^ N+l < e<5(l — 7). Then for n> N, 

(a) The polynomials P n (t) = Y^j=objt J have a unique root 7„ in [0, .6491] 

(b) Any 7 E [ji ow , 'Jhigh] which satisfies \Pn(l)\ < (liow) n also satisfies (7 - 7] < C\j\ n , 



where C 



1 

,5(1-7) • 



This theorem applies to beta encoders, corresponding to implementing the recursion 
([2]) with flaky quantizer Q u {u) defined by Q, and with 7 = known a priori to be 
contained in an interval [7io«> > Ihigh 

]. If /3 > 1.54059, or jhigh — -6491, then we can recover 
7 = f3~ 1 from either a bitstream corresponding to 0, or a pair of bitstreams (x, —x), using 
Theorem (4.1 ). Of course we should not consider only the scheme Q, but rather a revised 



scheme which accounts for integrator leak on the (single) integrator used in the beta- 
encoder implementation (see Figure 1). The revised beta-encoding scheme, with slightly 
different initial conditions, becomes 

Ul = x 
bi = Q(ux) 
for j > 1 : Uj + i = Xf3(uj — bj) 

b j+1 = Q(u j+1 ) (20) 

where A is an unknown parameter in [.9,1]. As long as (3 = Xj3 > 1, we still have that 
\x — Si=o^+l.7 i | < Cpj3~ N where = J-j-; furthermore, we can use Theorem (4.1) to 

recover /3 = A/3 in (20), although the specific values of A and f3 cannot be distinguished, 
just as the specific values of Ai and A2 in the expression for 7 could not be distinguished 
in GRE. 
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4.1 Remarks 



Figure 2 suggests that the first positive roots of the P n {t) serve as exponentially precise 
approximations to 7 for values of 7 greater than .6491; Figure 4 suggests that the first 
positive root of P n {t) will approximate values of 7 up to 7 = .75. Furthermore, these 



figures suggest that the constants C\ and C2 of (17) in the exponential convergence of 
these roots to 7 can be made much sharper, even for larger values of 7. This should 
not be surprising, considering that nowhere in the analysis of the previous section did 
we exploit the specific structure of beta-expansions obtained via the particular recursions 
(20) and ([7]), such as the fact that such sequences (b n ) cannot contain infinite strings of 



successive l's or -l's. It is precisely power series with such infinite strings that are the 



"extremal cases" which force the bound of .6491 in Theorem (3.1 ). It is difficult to provide 



more refined estimates for the constants C\ and C2 of (17) in general, but in the idealized 
setting where beta-expansions of are available via the ideal GRE scheme (J7| without 
leakage, or via the beta-encoding Q with /3 = <j>, the beta-expansions of have a very 
special structure: 

Theorem 4.2. Consider the ideal GRE recursion ([7]) with input uq = u\ = and 
Q(u,v) = Q u a (u,v), or the beta-encoder recursion ^ with f3 = (j), uq = 0, and Q(u) = 
Q u (u). As long as a > v in Q, or v < 1 in then for each j G {0, 1, ...}, bsj is equal 
to -1 or +1, and b 3j+ i = 63,4-2 = 

The proof is straightforward, and we omit the details. 



Theorem (4.2) can be used to prove directly that 7 must be the first positive root of 
the polynomials = &i + S?=i bj+itf , when 7 = </> _1 in either (JjJ or @. Indeed, 

-P/v(t) can be factored as follows: 

N 

PnW = ^(&3j+l - b 3j+1 t - b 3j+1 t 2 )t 3 i 
3=0 

N 

= (1 - i - i 2 ) ^>;+ii 3j 

= (1 - t - t 2 )R N (t) (21) 
where i?jv(t) is a polynomial with random coefficients of the form ^2jLq ±i 3 -' . 

P/v(t) = (1 — t — t 2 )i?Ar(t) clearly has a root at t = (f)^ 1 , and this root must be the 
only root of P/v(i) on [0,^ _1 ), since on this interval Rjsi{t) is bounded away from 
by |-Rjv(i)l > 1 — > 1 — Apa ~ -691. We can also obtain a lower bound on 

I^OT 1 )! = "(I + 2c/ ) - 1 )R N (0- 1 ) by IP^OT 1 )! > |1 + 20- 1 ||l - i^l > 1-545. Note 
that this bound holds also for the infinite sum F(t) = b\ + Yl'jLi bj+it 3 , and that the 
bound of |-F'(7)| > 1.545 is much sharper than the bound on |-F'(7)| given by Theorem 



(3.1); e.g., (5.63 = -07, and <5.649i = .00008 (see [I]). Similar bounds on the derivatives 



|-Pjv(7)| and |i^(7)| corresponding to beta-expansions of in a base 7 close to </> 1 should 



hold, leading to sharper estimates on the constants C\ and C2 of (17) in the case where 
beta-expansions of are available. 
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5 Conclusion 



In this paper, we have shown that golden ratio encoders are robust with respect to leakage 
in the delay elements of their circuit implementation. Although such leakage may change 
the base 7 in the reconstruction formula \y — Ylj=i ~ 0(j N ) , we have shown that 
exponentially precise approximations 7 to 7 can be obtained from the bitstreams of a pair 
(x,— x), and such approximations 7 are sufficient to reconstruct subsequent input y by 

|y-Ef=i 6^1-0(7^)- 

Our method can be extended to recover the base /3 in general beta-encoders, as long 
as (5 is known a priori to be sufficiently large; e.g. (3 > 1.54. Our method is similar to 
the method proposed in [3] for recovering (3 in beta-encoders when {0, l}-quantizers are 
used, except that our method does not require a fixed reference level, which is difficult to 
measure with high precision in practice. 




5 10 15 20 25 30 5 10 15 20 25 30 



(a) (b) 

Fi gure 2: (a) iV versus I7— 7/v|, where 7/v is obtained via a 10 step Newton Method approximation 
to the first root of the polynomial Pjv starting from xq = .618. For each N, I7 — jn\ is the worst- 
case error among 100 experiments corresponding to different (xj, —Xj) pairs and different 7^ chosen 
randomly from [.618, .7]. (b) iV versus \x — xn\, where xn is reconstructed from 7/v using the first 
N bits. The quantizer used in this experiment is Q v a (u, v) with v = .3 and a chosen randomly in 
the interval [1.7,2]. 



6 Appendix: Proof of Theorem (2.2) 



In this section we prove Theorem (2.2), which provides a range within which the "flak- 



iness" parameter v and the "amplifier" parameter a in the quantizer Q^(u,v) can vary 
from iteration to iteration, without changing the fact that the sequences (u n ) produced 



by the scheme (12) will be bounded. The derived range for v and a is independent of the 
specific values of the parameters (Ai,A2) € [.9, l] 2 in (12). 



The techniques of this section are borrowed in large part from those used in [5] to prove 
a similar result for the ideal GRE scheme ([7]); i.e., taking Ai = A2 = 1 in (23 ). As is done 



in [5 a , we first observe that the recursion formula (12) corresponding to the leaky GRE is 
equivalent to the following piecewise afiine discrete dynamical system on M 2 : 
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(a) (b) (c) 

Figure 3: The graphs of (a) Pg, (b) P\e, and (c) P32 for several pairs (xj, —Xj) corresponding to 
7 = .64575. These polynomials can have several roots on the unit interval, but the first of these 
roots must become exponentially close to 7 as N increases. The quantizer used is Q^(u,v) with 
parameter values v = .3, and a — 2. 




(a) (b) (c) 

Figure 4: The graphs of (a) P$, (b) Pi6, and (c) P32 for several pairs (xj, — Xj) corresponding to 7 = 
.75. These graphs suggest that the first positive root of these polynomials becomes exponentially 
close to 7 as TV increases, although we do not have a proof of this result for values of 7 greater 
than approximately .65. The quantizer used is Q^(u,v) with parameter values v = .3, and a — 2. 
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Figure 5: The rectangle R = Af B-2Cf Di (bold outline) is a positively invariant set for the map 
T£. In fact, T£(R) is contained in the smaller rectangle ABCD (dashed outline), illustrating that 
the revised GRE scheme (12 1 is robust with respect to small additive errors. Here, <&i and $2 are 
the eigenvectors of the matrix in A in (26) below. In this figure, Ai = A2 = .96, and \x is taken to 
be .0625. The length parameters l,r,h, and d will be discussed later. 



where 



Un+l 
Un+2 



U„ 
Un+l 



(22) 



rpU . 



We will construct a class of subsets R 



u 




1 




u 


~ Qa(u,v) 


' " 


V 




A1A2 Ai 




V 


1 



(23) 

R(n) = i?(Ai,A 2 ,M) of R 2 for which T»(R) + 
B u (0) C R, where B a (0) is the disk of radius \i centered at the origin; i.e. B^O) = 
{(u,v) : u 2 + v 2 < /U 2 }. Note that these sets are not only positively invariant sets 
of the map T£, but have the additional property that if (uq,vq) G R((j), the image 
(u n+ i, v n+ \) = T^(u n ,v n ) may be perturbed at any time within a radius of // (for example, 
by additive noise), and the the resulting sequence (u n ) n=0 will still remain bounded within 
R{jj) for all time n. 



We refer the reader to Figure 5 as we detail the construction of the sets R([i). Rect- 
angles A\B\C\D\ and A2B2C2D2 in Figure 5 are designed so that their respective images 
under the affine maps T\ and T2 (defined below) are both equal to the dashed rectangle 
ABCD. 



Ti : 
T 2 : 



u 




1 




u 




' " 


V 




A1A2 Ai 




V 




1 



u 

V 





AiA 2 



1 

Ai 



u 

V 



+ 



(24) 
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That is, rectangles A x BxC x D x and A2B2C2D2 will satsify Ti(A x B x CiDi) = ABCD = 
T 2 {A 2 B 2 C 2 D 2 ). More specifically, T X (A X ) = T 2 {A 2 ) = A, T X (B X ) = T 2 {B 2 ) = B, and so 

kL kh 

on. Since the rectangle = A^B^^ D\ is contained within the union of A X B X C\D X and 
A 2 B 2 C 2 D 2 , R is a positively invariant set for any map T(u, u ) satisfying 

r Ti(n,w), (u, v) £ R \ A 2 B 2 C 2 D 2 

T(u,v) = l T 2 (u,v), (u, v) £ R \ AiBiC\Di (25) 

( Ti(u,w) or T 2 (u,v), (u,v) e ff 

where -ff = j4ii?iCi-Di n A 2 B 2 C 2 D 2 . In particular, if the parameters a and z/ in the map 
are chosen such that the intersection of R and the strip F = {(it, u) : — f < u + av < u} 



is a subset of i/, then is of the form (25 ). Indeed, F is the region of the plane in which 
the quantizer Q^(u,v) operates in flaky mode. This geometric setup is clarified with a 
figure, which we provide in Figure 6. 




Figure 6: The flaky quantizer Q^(u, v ) outputs either —1 or +1 when the input (u, v) belongs to 
the strip F = {(u,v) : —v < u + av < u} (shaded above). Here \i — .0625, a = 2, and v = .15. 
For these parameter values, the intersection of F and R is a subset of A\B\C\D\ D A2B2C2D2, 



meaning that the scheme (12 1 produces bounded sequences (u n ), when implemented with this 



quantizer, and with input (uq,Ui) G R. 

It remains to verify the existence of at least one solution to the setup in Figure 5 
for each (Xi,X 2 ) £ [.9, l] 2 . This can be done, following the lead of [5], in terms of the 
parameters defined in Figure 5. 



Note that the matrix 



A 





A1A2 



1 

Ai 



in (|24|) has as eigenvalues ei 
A 2 



Ai+VAf+4A 1 A 2 



and — e 2 



when Ai = A2 = 1, we have that ei = <f> w 1.618 and e 2 
and —e 2 have respective normalized eigenvectors 



1 

ei 



and $2 



-1 



(26) 



-Ai+^/Af-4AiA 2 



) . In particular, 
.618. The eigenvalues ei 



1 

(2 
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where si = Wl + ef and S2 = \/l + ef- It follows that the affine map Ti acts as an 
expansion by a factor of ei along $1 and a reflection followed by contraction by a factor 
of e 2 along $2; followed by a vertical translation of +1. T 2 is the same as T\ except with a 
vertical translation of -1 instead of +1. After some straightforward algebraic calculations, 
the mapping relations described above imply that the parameters in Figure 5 are given by 

2/i 2si 
h = h 



1-ei ei(ei - l)(ei + e 2 ) 
d = -^- + Sl(2 " £l) 



1-ei ei(ei - l)(ei + e 2 ) 
j . /f , s s(2 - e 2 ) 



l-e 2 e 2 (l - e 2 )(ei + e 2 ) 

" + * 5 . (27) 



1 - e 2 (1 - e 2 )(ei + e 2 ) 

The existence of the overlapping region H = A\B\C\D\ D A2B2C2D2 is equivalent to the 
condition d > which in turn is equivalent to the condition /i < ^j^p^A • This expression 
is minimized over the range (Ai, A 2 ) € V = [.9, l] 2 when Ai = A 2 = 1, in which case this 
constraint simplifies to « < ~ -2008. 

We are interested primarily in quantizing numbers in [—1,1], so that this is the range 
of interest for the input u\ = x; for no we simply take 0. The set {0} x [— 1,1] c 
H(Ai x 2 )ev ^(^l) A*) is equivalent to the condition that the line v = mu + b which passes 
through the points B and C in Figure 5 have y- intercept b > 1 for all (Ai, A 2 ) G V. This 
line has slope m = — e 2 , and passes through the point C, so that its y-intercept is given 
by b = s^ 1 (ei + e 2 )(/t — d). Rearranging terms, this implies that b > 1 if and only if 

ei + e 2 V ^ 

The right hand side of the above inequality is bounded below over the range (Ai, A 2 ) G V 
by V 1 +(- 9 ^ 2 ( 2 -^ ) ^ 301) so that in deed {0} x [-1, 1] c R(Xi, A 2 , u) is satisfied indepen- 
dent of (Ai, A 2 ) G T/, for all admissable \i G [0, .2008]. 

In practice, the "flaky" parameter v of the quantizer Q v a is not known exactly; we will 
however be given a tolerance 5 for which it is known that \i>\ < 5. It follows that for each 
5, we would like a range (a m m, a m aa;) f° r the "amplifier" a such that the GRE scheme 



(12), implemented with quantizer Q v a , produces bounded sequences (u n ) for all values of 
(Ai,A 2 ) G V. For a particular choice of (Ai,A 2 ), it is not hard to derive an admissable 
range for a in terms of the eigenvalues e\ and 62 for fixed 5. We will take u = in the fol- 
lowing analysis for the sake of simplicity. First of all, the tolerance 5 must be admissable, 
i.e. the coordinate (6,0) should lie within the region H. Indeed, for v G [0,(5], a can vary 
within a small neighborhood of ^ , as long as the line with slope — ^ which passes through 
(0, 6) remains bounded above in the region H by the line passing through B\ and C\. In 
other words, the admissable range for a is obtained from the constraint m\ < — ^ < m 2 , 
where mi is the slope of the line through the points {(6,0), C±} and 7712 is the slope of 
the line through the points {(6, 0), B^j in Figure 6. Rewriting C\ and B^ in terms of the 



eigenvalues ei and 62 via the relations (27), we have that 
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L(e 1 ,e 2 ) <a< ?7( e l> e l + e 2) 



(29) 



where 



and 



L(x,y) = N(x,y)/D(x,y) 

x(x - 1) - (2 - x)(l - y) + &c(x - l)(x + y)(l - y) 
x((2-x)(l-y) + y(x-l)) 



tt( \ ^ + ^y -2y - 8xy(x - 1)(1 - y + x) 

U(x, y) = -j—^ . (30) 



It is not hard to show that for any fixed y G [.90 + .90 _1 , <j) + (f)~ 1 } ps [2.0124, 2.236], the 
function U(x,y) increases as a function of x over the range x G [.90,0] ~ [1.456,1.618], 
as long as v < .75. It follows that the minimum of U{x,y) over the rectangle S = 
[1.456,1.618] x [2.0124,2.236] occurs along the edge corresponding to x = 1.456. But 
£7(1.456, y) decreases as a function of y over the interval [2.0124,2.236], so that the min- 
imum of U(x,y) over the entire rectangle S occurs at (x,y) = (1.456,2.236), giving the 
following uniform lower bound on o. max for (Ai,A 2 ) G V: 

amax > [7(1.456,2.236) 

w 2.281 - .9525. (31) 

We proceed in the same fashion to derive a uniform upper bound on a m i n over (Ai, A2) G 
V, except that we analyze the numerator and the denominator in the expression for 
L(x, y) separately. The numerator N(x, y) increases as a function of x over the range 
x e [1.456,1.618] for fixed y G [.90 -1 ,<£ -1 ] ps [.5562, .618], so that N(x,y) achieves its 
maximum over the rectangle T = [1.456, 1.618] x [.5562, .618] along the edge corresponding 
to x = 1.618. N (1.618, y) increases as a function of y over the range y G [.5562, .618], so 
that N(x,y) achieves its maximum over T at (x,y) = (1.618, .618). A similar analysis 
shows that the denominator D(x,y) is minimized over T at (x,y) = (1.456, .618). It 
follows that for (Ai,A 2 ) G [.9,1] 2 , 

(Xmin is bounded above by 

amin < iV(1.618, .618)/£>(1.456, .618) 

« 1.198(1 + 5) (32) 

Finally, we note that 5 is admissable for all (Ai, A 2 ) if and only if L(e±, e 2 ) < U(e\,e\ + e 2 ) 
for this value of 8 and for all coresponding e\ and e 2 . Thus an admissable range for 8 is 



obtained by equating the lower bound of 2.281 — .9525 in (31) and the upper bound of 



1.198(1 + 8) in (32); namely, 8 G [0, .5037]. 
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